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Abstract 

Stochastic waveforms are constructed whose expected autocorrelation 
can be made arbitrarily small outside the origin. These waveforms are 
unimodular and complex-valued. Waveforms with such spike like auto- 
correlation are desirable in waveform design and are particularly useful in 
areas of radar and communications. Both discrete and continuous wave- 
forms with low expected autocorrelation are constructed. Further, in the 
discrete case, frames for are constructed from these waveforms and the 
frame properties of such frames are studied. 
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1 Introduction 
1.1 Motivation 

Designing unimodular waveforms with an impulse-like autocorrelation is cen- 
tral in the general area of waveform design, and it is particularly relevant in 
several applications in the areas of radar and communications. In the former, 
the waveforms can play a role in effective target recognition, e.g., [1], [2], [3], 
[4], [5], [6], [7], [8]; and in the latter they are used to address synchronization 
issues in cellular (phone) access technologies, especially code division multiple 
access (CDMA), e.g., [9], [10]. The radar and communications methods com- 
bine in recent advanced multifunction RF systems (AMRFS). In radar there 
are two main reasons that the waveforms should be unimodular, that is, have 
constant amplitude. First, a transmitter can operate at peak power if the signal 
has constant peak amplitude - the system does not have to deal with the sur- 
prise of greater than expected amplitudes. Second, amplitude variations during 
transmission due to additive noise can be theoretically eliminated. The zero 
autocorrelation property ensures minimum interference between signals sharing 
the same channel. 
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Constructing unimodular waveforms with zero autocorrelation can be re- 
lated to fundamental questions in harmonic analysis as follows. Let R be the 
real numbers, Z the integers, C the complex numbers, and set T = M/Z. The 
aperiodic autocorrelation Ax : Z — )• C of a waveform X : Z — )• C is defined as 

1 ^ 

VfceZ, Ax[k] = lirn^^^ X[k + m]X[m]. (1) 

m=-N 

A general problem is to characterize the family of positive bounded Radon mea- 
sures F, whose inverse Fourier transforms arc the autocorrelations of bounded 
waveforms X. A special case is when F = 1 on T and X is unimodular on Z. This 
is the same as when the autocorrelation of X vanishes except at 0, where it takes 
the value 1. In this case X is said to have perfect autocorrelation. An extensive 
discussion on the construction of different classes of deterministic waveforms 
with perfect autocorrelation can be found in [11]. Instead of aperiodic wave- 
forms that arc defined on Z, in some applications, it might be useful to construct 
periodic waveforms with similar vanishing properties of the autocorrelation func- 
tion. Let n > 1 be an integer and Z„ be the finite group {0, 1, . . . , n — 1} with 
addition modulo n. The periodic autocorrelation Ax : Z„ ^ C of a waveform 
X : Z„ — >• C is defined as 

^ n—l 

VA; = 0,l,...,n-l, Ax[k] = - ^ X[m + k]X[m]. (2) 

It is said that X : Z„ — > C is a constant amplitude zero autocorrelation 
(CAZAC) waveform if each \X[k]\ = 1 and 

n—l 

Vfc = l,...,n-1, Ax[k] = -Y,X[^ + k]X\m\ = 0. 

^ m=0 

The literature on C AZ AGs is overwhelming. A good reference on this topic is [3] , 
among many others. Literature on the general area of waveform design include 
[12], [13], [14]. Comparison between periodic and aperiodic autocorrelation can 
be found in [15]. 

Here the focus is on the construction of stochastic aperiodic waveforms. 
Henceforth, the reference to waveforms shall imply aperiodic waveforms unless 
stated otherwise. These waveforms are stochastic in nature and are constructed 
from certain random variables. Due to the stochastic nature of the construction, 
the expected value of the corresponding autocorrelation function is analyzed. It 
is desired that everywhere away from zero, the expectation of the autocorrela- 
tion can be made arbitrarily small. Such waveforms will be said to have almost 
perfect autocorrelation and will be called zero autocorrelation stochastic wave- 
forms. First discrete waveforms, AT : Z C, arc constructed such that X has 
almost perfect autocorrelation and for all n G Z, |A[n]| = 1. This approach is 
extended to the construction of continuous waveforms, a; : K — >■ C, with similar 
spike like behavior of the expected autocorrelation and \x{t)\ = 1 for all t gR. 
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Thus, these waveforms are unimodular. The stochastic and non-rcpetitivc na- 
ture of these waveforms means that they cannot be easily intercepted or detected 
by an adversary. Previous work on the use of stochastic waveforms in radar can 
be found in [16], [17], [18] where the waveforms arc only real- valued and not 
unimodular. In comparison, the waveforms constructed here are complex val- 
ued and unimodular. In addition, frame properties of frames constructed from 
these stochastic waveforms are discussed. This is motivated by the fact that 
frames have become a standard tool in signal processing. Previously, a mathe- 
matical characterization of CAZACs in terms of finite unit-normed tight frames 
(FUNTFs) has been done in [2]. 

1.2 Notation and mathematical background 

Let X be a random variable with probability density function /. Assuming X 
to be absolutely continuous, the expectation of X, denoted by E{X), is 

E{X) = / xf{x) dx. 
Js. 

The Gaussian random variable has probability density function given by f{x) = 
— ^e^^'^^s^' . The mean or expectation of this random variable is u and the 

variance, V{X), is a^. In this case it is also said that X follows a normal dis- 
tribution and is written as X ^ N{^,a'^). The characteristic function of X 
at t, i5(e**'^), is denoted by (j>xit)- For further properties of expectation and 
characteristic function of a random variable the reader is referred to [19]. 

Let H he a Hilbert space and let V = {vk, k € JC}. whcirc; /C is some index 
set, be a collection of vectors in H. Then V is said to be a frame for H if there 
exist constants A and B, < A < B < oo, such that for any v & % 

A\\vr<Y,\{v,Vk)?<B\\vr. 
keic 

The constants A and B are called the frame bounds. Thus a frame can be 
thought of as a redundant basis. In fact, for a finite dimensional vector space, 
a frame is the same as a spanning set. li A = B, the frame is said to be tight. 
Orthonormal bases are special cases of tight frames and for these, A = B = 1. 

If is a frame for H then the map F -.H ^ (-2{^) given by F{v) = {{v, Vk) : 
/e e /C} is called the analysis operator. The synthesis operator is the adjoint 
map F* : ^2(/C) H, given by 

P*{{0'k}) = "^akVk- 

keK. 

The frame operator : "H — > "H is given by = F* F. For a tight frame, the 
frame operator is just a constant multiple of the identity, i.e., = AX, where 
I is the identity map. Every v gH can be represented as 

v = ^{v,J^~^Vk)Vk = ^{v,Vk)J^~'^Vk. 

keic keic 
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Here {T~^Vk} is also a frame and is called the dual frame. For a tight frame, 
J^~^ is just ^J^. Tight frames are thus highly desirable since they offer a com- 
putationally simple reconstruction formula that does not involve inverting the 
frame operator. The minimum and maximum eigenvalues of T are the optimal 
lower and upper frame bounds respectively [20]. Thus, for a tight frame all the 
eigenvalues of the frame operator are equal to each other. For the general theory 
on frames one can refer to [20], [21]. 

1.3 Outline 

The construction of discrete unimodular stochastic waveforms, X : Z — >• C, 
with almost perfect autocorrelation is done in Section 2. This is first done 
with the Gaussian random variable and then generalized to other random vari- 
ables. The variance of the autocorrelation is also estimated. The section also 
addresses the construction of stochastic waveforms in higher dimensions, i.e., 
construction oi v : Z that have almost perfect autocorrelation and are 

unit-normed, considering the usual norm in C''. In Section 3 the construction of 
unimodular continuous waveforms with almost perfect autocorrelation is done 
using Brownian motion. 

As mentioned in Section 1.2, frames are now a standard tool in signal process- 
ing due to their effectiveness in robust signal transmission and reconstruction. 
In Section 4, frames in C'' (d > 2) are constructed from the discrete waveforms 
of Section 2 and the nature of these frames is analyzed. In particular, the maxi- 
mum and minimum eigenvalues of the frame operator are estimated. This helps 
one to understand how close these frames are to being tight. Besides, it follows, 
from the eigenvalue estimates, that the matrix of the analysis operator, for 
such frames, can be used as a sensing matrix in compressed sensing. 

2 Construction of discrete stochastic waveforms 

In this section discrete unimodular waveforms, X : Z — >■ C, are constructed from 
random variables such that the expectation of the autocorrelation can be made 
arbitrarily small everywhere except at the origin. First, such a construction is 
done using the Gaussian random variable. Next, a general characterization of 
all random variables that can be used for the purpose is given. 

2.1 Construction from Gaussian random variables 

Let {Y^l^gz be independent identically distributed (i.i.d.) random variables 
following a Gaussian or normal distribution with mean and variance cr^, i.e., 
Yi ~ N{0, 0-2). Define X : Z C by 

Vn e Z, X[n] = ^"=-" (3) 
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where i is -/—I. Thus, for each n, l-'^fn]] = 1 and X is unimodular. The 
autocorrelation of X at A; G Z is 

where the hmit is in the sense of probabihty. Theorem 2.1 shows that the wave- 
form given by (3) has autocorrelation whose expectation can be made arbitrarily 
small for all integers fc ^ 0. 

Theorem 2.1. Given e > 0, the waveform X : Z ^ C defined in (3) has 
autocorrelation Ax such that 

Proof, (i) When fc = 0, 



-N 



and so E{AxW) = 1. 

(ii) Let fc > 0. One would like to calculate 



\ n=-N J 

Let gN{X) = ^Y.n=-NX[n + k]X\^]. Then \gN{X)\ < 1. Let h{X) = 1. 
Then for each TV, \gNiX)\ < h{X) and E[h{X)] = 1. Thus, by the Dominated 
Convergence Theorem [19], which justifies the interchange of limit and integra- 
tion below, one obtains 

E{Ax[k]) = sflun^^^ f; X[n + k]X[^]\ 

\ ^ n=-N J 

= lim — — - V E{X{n + k\X[n\) 



lim — -— — 

iV-)-oo IN + 

lim 



n=-N 

N 



- J2 E{e-^^l 



n=-N 

N 



n=-N 

N 



= lim ^ ^ ^ 
N^oo 2N + 



1 w— , r / n ■ \n2fc r /o- \n2fc 



n=-N 
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where the last hnc uses the fact that the Ygs arc i.i.d. random variables. Here 
(^) is the characteristic function at ^ of Yi which is the same as that for 
any other Yg due to their identical distribution. The characteristic function at 

— of a Gaussian random variable with mean and variance is e~ 
Thus 

E{A^[k]) = [e-4m'f = e-'^^'i'?r. 
(iii) When > 0, a similar calculation for fc] gives 

Together, this shows that given e and any fc 7^ 0, 

£;(^x[fc]) = e-l'=l<^'(^)' 

which indicates that the (;xpcctation of the autocorrelation at any integer fc 7^ 
can be made arbitrarily small depending on the choice of e. □ 

As shown in Theorem 2.1 the expectation of the autocorrelation can be made 
arbitrarily small but this is not useful unless one can estimate the variance of 
the autocorrelation. Denoting the variance of Ax [A;] by one has 

V{Ax[k]) - EQAxM') - \EiAxm^ = E{\Ax[k]\^) - e'^l'^l'^'^")'. 

First consider A: > 0; 

/ 1 ^ \ / 1 M 

lAxM"^ = lim — — - y X\n + k]X\n]] lim — — - V X\m + k]X\r. 

\ n=-N / \ m=-M 



1 



1 



N 



M 



= lim lim 

JV->oo M^cx> {2N + 1) (2M + 1) 

= lim lim 



JV->oo M->oo {2N + 1) (2M + 1) 



J2 J2 X[n + k]X[n] X[m + k]X[m] 

2= — N m=—M 
N M 

E E 



-N 1 



By applying the Lebesgue Dominated Convergence Theorem one can bring the 
expectation inside the double sum to get 



E{\Ax[k]\^) = lim lim 



N^oo M^oo {2N + 1) (2M + 1) 



1 ^ N M 

+ 1) (2M+1) E E 



n=—N m=—M 

E (^e^i{T.U Y-n-,+^U Yn+,-Y,U y-^-i-^U ^-+^)) 
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The sum 

k k k k 

^-n-j + ~ ^-m-j ~ ^m+j (4) 

i=i j=i i=i i=i 

may have cancelations among terms involving n with terms involving m. Sup- 
pose that for a fixed n and m, there are km,n indices that cancel in each of 
the 4 sums in (4). Due to symmetry, the same number i.e., kmm of terms will 
cancel in each sum. Depending on n and m, kmn lies between and fc, i.e., 
< kmn < k. For the sake of making the notation less cumbersome, kmn will 
from now on be written as k. When m~n, k~k.lim, >n + kovn>m + k 
then A: = 0. Each sum in (4) has k terms and k of these get cancelled leaving 
[k — k) terms. One can re-index the variables in (4) and write it as 



k k k k 

' VL , .• = + . . . + y. 

4()c-fe) 



i=i i=i i=i 



where the sign depends on whether m is less than or greater than n. Thus 



^' Jl jv^oo M^-oo (27V + 1) (2M + 1) 

Due to the independence of the Y^s, this means 

= lim lim = ^ V V SyA ±—\ 



N M 



N M 



4(fc-fe) 



n=—N m=—M 
1 -I AT M 

= lim hm , , V V e-'rm'^i^-'^) 

AT^-ooM^-oo (27V + 1) (2M + 1) ^ ^ 

^ ' ' ^ ' ' n=-N m^-M 

The minimum is attained for fc = and the maximum at k = k. Thus 

1 . N M 

E{\Ax[k]\^) < lim lim E E 1 = 1 and 

^1 Jl ^ - JV^-ooM^-oo (2iV-|- 1) (2M+ 1) 

^ ' V I / n=-N m=-M 

1 1 ^ ^ 2 •> 2 

i;(|A;,[fc]n > lim lim E E e-^(^)"4fc = e-V(¥) 



This gives 

0<n^x[fc])<l-e-2'=-'(^r. 
A similar calculation can be done for fc < 0. Thus for fc ^ 0, 

0<y(Ax[fc])<l-e-2|^l-'("r. 
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2.2 Generalizing the construction to other random Vciri- 

ables 



So far the construction of discrete unimodular zero autocorrelation stochastic 
waveforms has been based on Gaussian random variables. This construction 
can be generalized to many other random variables. The imimodularity of the 
waveforms is not affected by using a different random variable. The following 
theorem characterizes the class of random variables that can be used to get the 
desired autocorrelation. 

Theorem 2.2. Let {Y^j^gz be a sequence of i.i.d. random variables with char- 
acteristie function (py- Suppose that the probability density function of the Yf S is 
even and that 4>Y{t) goes to as t goes to infinity. Then, given e, the waveform 
X -.Z ^ C given by 



X[n]- 

has almost perfect autocorrelation. 



e « 



Ye 



Proof. Since the density fimction of each Yi is even this means that the char- 
acteristic function is real valued [19]. Following the calculation in the proof of 
Theorem 2.1, the expected autocorrelation of X for A; 7^ is 



E{Ax[k]) 



2\k\ 



and this goes to zero with e by the hypothesis. 



□ 



Example 2.3. Suppose the l^s follow a bilateral distribution that has density 
e~l^l with X e (—00,00) and characteristic function = Then for 

2|fe| 



and this can be made arbitrarily small with e. 

In the same way as was done in the Gaussian case, for A; > 0, 



E(\A 



Thus 



< 



lim lim ^ ^- 

jv^-oo M^oo (27V + 1) (2M H 

1 and 



1 



1 + 



4k 



1) 



N M 

E E 

i=—N m=—M 



4|fe| 



Q<V{Ax[k])<l- 



1 + 



27r 
± — 



4(fe-fe) 
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Example 2.4. Suppose that the Y^s follow the Cauchy distribution with density 
function ^p^j^^- Note that, disregarding the constant n, this is the character- 
istic function of the random variable considered in Example 2.3. The charac- 
teristic function of the Ygs is now e~l*l, the same as the distribution function in 
Example 2.3. For fc ^ 0, 

47r|fc| 

= e 

which can be made arbitrarily small with e. Also, 

0<V{Ax[k]) <l-e — —. 



E{Ax[k]) 



2.3 Higher dimensional case 



Here one is interested in constructing waveforms t; : Z — ^ C*, d > 2. It is desired 
that V has unit norm and the expectation of its autocorrelation can be made 
arbitrarily small. One way to construct v is based on the construction of the 
one dimensional example given in Section 2.1. This is motivated by the higher 
dimensional construction in the deterministic case [2]. As before, {Y^j^gz is 
a sequence of i.i.d. Gaussian random variables with mean zero and variance 



a . Next, one defines X[n] = e « 



defined as 



Vm e Z, v[m] = —j= 



1 



The waveform v 

X[m] 
X[m+l] 

X[m + d-l] 



is then 



(5) 



In this case, the autocorrelation is given by 



JV 



Ay\k] = lim — — — - 



^ {v[n + k],v[n]) 



(6) 



n=-N 



where (., .) is the usual inner product in C*. The length or norm of any v[m] is 
thus given by 

||v[m]||^ = (t;[m], u[m]). 

From (5), 

II^Mf = ^Ex[m + n]X[m + n] = ^ = l. 



ri=0 



Thus the v[m\s are unit-normed. The following Theorem 2.5 shows that the 
expected autocorrelation of v can be made arbitrarily small everywhere except 
at the origin. 
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Theorem 2.5. Given e > 0, the waveform u : Z — >■ defined in (5) has 
autocorrelation Ay such that 



1 i/ fc = 



E{Ay[k])-- 
Proof. As defined in (6), 

1 ^ 

^"[^^ = J^oo 27VTT ^ + ''^'^I"^^- 



When k = 0, 



n=-N 



N 



n=-N 



Thus, 

For fc ^ 0, due to (5), 



E{AM) = 1. 



{v[n + k],v[n]) 



1 



X[n + k] 
X[n + k + l] 



X[n] 
X[n + 1] 

X\n + d- 11 



X[n + fc + d- 1] 
= i {x[n + + X[n + A: + l]X[n + 1] + . . . + 

X[n + A; + d- l]X[n + rf - 1]). 



Consider > 0. 



E{Ay[k]) = ^Imi^^^ E i?((^[n + fc],t;[n])) 

= lim — — - V E\ X\n + k + m]X[n + 

n=-N \ m=0 



N^oo 2N + 



n— — N 7n—0 
^ N d~l 



^i(F_(„+„+fc) + ...+F_(„+„+i)+Y„+„+i+...+Y„+„+fc)^ 



n— — N m—0 
N , d-l 



N^oo 2N + 



1 y l^yEie'^^^^r = e-^'W. 



n— — N m—0 
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Similarly, for A; < one gets 

□ 

Thus the waveform v as defined in this section is unit-normed and has autocor- 
relation that can be made arbitrarily small. 

Remark 2.6. As in the one dimensional construction, it is easy to see that 
here too the construction can be done with random variables other than the 
Gaussian. In fact, all random variables that can be used in the one dimensional 
case, i.e., ones satisfying the properties of Theorem 2.2, can also be used for the 
higher dimensional construction. 



2.4 Remcirk on the periodic CEise 

It can be shown that the periodic case; follows the same nature as the aperiodic 
case. The sequence X : Z„ — >■ C is defined in the same way as in Section 2.1, 
i.e., 

Vme {0,l,...,n- 1}, X[m] = e^'^'^=-m^^ 
where N{0,a'^). Following the definition given in (2), when k = 0, 



^ n — 1 



m=0 



When k ^0, the expectation of the autocorrelation is 



n— 1 

E{Ax[k]) = -J2E{X[m + k]XH) 



m=0 



For fc > 0, 

E{Ax[k]) 



^ n— i 



TO=0 

n-1 



^ m=0 
^ n— 1 

TO=0 





2k 




)] 













where one uses the fact that the Y^s are i.i.d.. A similar calculation for negative 
values of k suggests that the autocorrelation can be made arbitrarily small, 
depending on e, for all non-zero values of k. Also, as in the aperiodic case, this 
result can be obtained for random variables other than the Gaussian. 
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3 Construction of continuous stochastic wave- 
forms 

In this section continuous waveforms with almost perfect autocorrelation are 

constriictcd from a one dimensional Brownian motion. 

For a continuous waveform a: : R — C, the autocorrelation A^: : M C can 
be defined as 

1 

Ms) = lim — / x{t + s)xit) dt. (7) 

Let {W{t); i > 0} be a one dimensional Brownian motion. Then W{t) satisfies 

(i) WiO) = 

(ii) W(t + s) - W{s) ~ 7V(0, (jH) 

(iii) < ti < ■ ■ ■ < tk, W{ti+i) — W{ti) are independent random variables. 

Theorem 3.1. Let W{t) be the one dimensional Brownian motion and e > 

be given. Define a; : M — > C by 

x{t)=e^'^^'^ fort>0, 
and x{—t) = x{t). Then the autocorrelation of x, A^, satisfies 

f 1 ifs = 

E{A^{s)) = i -4uU2^Y -f un 

Proof. We would like to evaluate 

E{A,{s)) = E (^lirn^ ^ J^^ x{t + s)^ dt^ . 

Let s > and let 5t(s) = ^ /Tt ^{^ + s)x{t) dt. 

1 

E{9t) = ^E{x{t + s)x{t)) dt 

= ^ j ^4>w{t+s)-w(t)(^^ =e-^''y^) <oo. 

Thus each qt is integrable and further \gT\ < 1- Let h{t) = 1; t e M. Then 
E{h) = 1. Therefore, by the Dominated Convergence Theorem, and properties 
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of Brownian motion and characteristic functions, one gets 



= lim — r E fe^ ^(^(*+«)-^(*))) dt 



= lim 



1 



''W(t+s)-Wit) 



27r 



T-)-00 2T J_rp 

which can be made arbitrarily small based on e. Similarly, 



E{A{-s)) = lim^ — I (l)w(t)-w{t-s) { - — 



27r 



= e 2 



□ 



4 Connection to frames 

Consider the mapping v : Z ^ given by 



v{k) = 



Vd 



( m \ 

X[k + l\ 



\ X[k + d-l\ ) 



(8) 



where X[k\ = e^' '^'^'=-'''^' , as defined in Section 2.1. 

Let M > d and consider the set V = {v{l),v{2), . . . ,v{M)} of M unit 
vectors in C^. The matrix 



Vd 



m x[2] 

X[2] X[3] 



m 

X[d+1] 



X[M] X[M+1] ■■■ X[M + d-l] 



is the matrix of the analysis operator corresponding to V. The frame operator 

of y is J" = F*F, i.e., 



^=-d 



X[l] X[2] 
X[2] X[A\ 



X[M] 
X[M + l] 



X[d] X[d+1] ■■■ X[M + d-l\ 



X[l\ X[2] 



X[2] X[3] 



X[d\ 
X[d+l] 



X[M] X[M+l\ ■■■ X[M + d-l] 
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The entries of T are given by T-m.-m = ^ and for m 7^ n, m > n, 



J~ rr 



-(X[m] X[m + 1] ■■■ X[M + m-l] 



/ X[n] \ 
X[n+l] 

\ X[M + n-l\ I 



M 1 



M-l 



x[m + e\x[n + e\ 




'-i(Y.^.i + -+Y.n-t--L+Yn+t+1 + -+Yra+e) 



d \ M 



Note that since F is self-adjoint, Tm.n = ^n,m- It is desired that V emulates 
a tight frame, i.e, T is close to a constant times the identity, in this case, ^ 
times the identity. Alternatively, it is desirable that the eigenvalues of T are all 
close to each other and close to In this case, due to the stochastic nature of 
the frame operator, one studies the expectation of the eigenvalues of T. 



4.1 Frames in 

This section discusses the construction of sets of vectors in as given by (8). 

The frame properties of such sets are analyzed. In fact, it is shown that the 
expectation of the eigenvalues of the frame operator are close to each other, the 
closeness increasing with the size of the set. The bounds on the probability of 
deviation of the eigenvalues from the expected value is also derived. The related 
inequalities arise from an application of Theorem 4.1 [22] below. 

Theorem 4.1 (Azuma's Inequality [22]). Suppose that {Xj. : fc = 0, 1, 2, . . .} 
is a martingale and 

\Xk - Xk-i \ < Cfe, 
almost surely. Then for all positive integers N and all positive reals t, 

P{\Xn - Xo\ >t)< 2e(^^Sfc^) . 
Consider M > 3 vectors in C^, i.e., d = 2 in (8). Then w : Z and 

(9) 
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Considering the set V = {v{l), v{2), . . . , v{M)}, the frame operator of V is 



X[l] X[2] X[3] 
X[2] X[3] X[4] 



X[M] 
X[M + 1] 



X[l] X[2] 



X[2] X[3] 



m x[4] 



or, 



M 
2 



^i:Z=iX[m]X[m+l] 



ilY.Z=iX[m]X[m+l] 



X[M] X[M + l] 

■ (10) 



Theorem 4.2. (a) Consider the set V = {«(!), w(2), . . . , v{M)} C C^, M > 3, 
where the vectors v{n) are given by (9). The minimum eigenvalue, Amin 
and the maximum eigenvalue, Amax(.^), of the frame operator ofV satisfy 



^{l-6)< E{X^i4J^)) < i^(Ama.(J^)) < Y (1 + (5) 



(11) 



where 6 = ^JT^IT^^W. 

(b) The deviation of the minimum and maximum eigenvalue of T from their 
expected value is given, for all positive reals r, by 

P(|A,„i„(J-) - ^(A„,i„(J-))| > r) < 2e-^, 

P(|A„>ax(.F) - £;(Aniax(.F))| > r) < 2e-^. 

Proof, (a) The frame operator of V = {v{l),v{2), . . . ,v{M)} is given in (10). 
The eigenvalues of -^T are Ai = 1 — \a\ and A2 = 1 + \a\ where 

M 



m=l 



Let 



71 = X[1]X[2] = 6-^-2+^2)^ 

72 = X[2]X[3] = e-" 



■i(Y-_3+y3) 



so that 



7m 



X[M]X[M + 1] = e-^'^^-f^+D+'^c^+i)^ 
71 + 72 H 1- 7m 



a 



M 



Note that for m ^ n, 7^ and 7„ are independent and so E{jrnjn) — E{^rn)E{j^). 
Also, since the Y^s are i.i.d. and the characteristic function of the Y^s is sym- 
metric, 



V 1 < m < M, ^(7„) = i;(e'r*(i'-(".+i)+^".+i)) = {£;(e^*'*'i)}2 = e 



E{lm) 
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and therefore 

Thus 



Eijm'ln) = e 



E{\a\^) = £;(aa) = -^£((71 +72 + ---+7M)(7r + ^+---+W)) 
= li^^(l7i|' + l72|' + --- + |7M|'+ ^7m7^) 



M 
1 

M ' Af2 



1 1 _s 



1 M-l 9„2,2^^2 

= 1 e 

The above estimate on iJ(|a|^) impUes that 



E{\a\)<^E{\o^) = \^,+ 



M-l 



M M 

Since £;(Ai) = 1 - E{\a\) and E{\2) = 1 + (12) implies 



(12) 



M M 



1 M-l 



M M 



Noting that XiauA^) = Ainax(^) = '^^2, one finally gets, after setting 



" V M ^ M ' 



^ (1 - (5) < ii;(A^i„(^)) < E{\^,^{F)) <^{l + 5). 

(b) To prove (b) wc use the Doob martingale and Azuma's incquahty [22]. For 
n = 2, . . . , M+1, let Zn-i = Y-n+Yn- Here the Doob martingale is the sequence 
{Uo, Ui,..., Um-i} where 



and 



for fc = 1,...,M- 1, Uk = E 



1 



1 

M 



M 



M 



\Zi,Z2, • . . , ^fe 



Note that i7o = E{\a\) and J/m-i = Also, 

|C/fc-C/fc-i|<|C/fe| + |C/fc-i|<2. 
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So by Azuma's Inequality (see Theorem 4.1) 



PQUm-i - Uo\ >r)= P{\\a\ - E{\a\)\ > r) < 2e ^^ET^ = 2e-^ . 
Since |Ai - E{Xi)\ = IA2 - E{X2)\ = ||q;| - E{\a\)\, this means 

P(|Ai - i;(Ai)| > r) < 2e-^ 

and 

P(|A2-i;(A2)| >r)<2e-i^. 

Going back to the actual frame operator whose eigenvalues are ^Ai and 
^A2, the following estimates hold. 

P(|A„,ax(^)-i?(A„,ax(^))| >r) = P(^y ||a|-i^(|a|)| >r 

= p(^\\a\-E{\a\)\>^r^<2e-f^ 

and 

P(|A„,i„(7-)-£(A^in(7-))| >r) = p(^^\\a\-E{\a\)\>r^ 

= P(^||a|-ii;(|a|)| > ^r) < 2e-^. 

□ 

Corollary 4.3. The eigenvalues of the frame operator considered in Theorem 
4-2 satisfy, for all positive reals r, 

P {x^Un < Y (1 - ^) - r) < e-4'-^/«^^ 

P (A^a.(^) > f (1 + <5) + r) < 



Proof. Due to part (a) of Theorem 4.2 

M 

Amin(^) < - (5) - r- ^ A„,in(-F) < £(A„,n(^)) - r. 

This impUes, as a consequence of part (b) of Theorem 4.2, that 

P < y (1 - <5) - < P (A^in(7-) < E{\^^{T)) - r) < g-^-VsM^^ 
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In a similar way, from part (a) of Theorem 4.2 



M 

Amax(J") > y(1 +5)+r =^ A„jax(J") > E{X^^{J^)) + r 
which implies, as a consequence of part (b) of Theorem 4.2, that 

P (A^ax(J^) > y (1 + 5) + < P > E{X^^{J^)) +r)< e-^-'m\ 

□ 

Remark 4.4. In Theorem 4.2, as M tends to infinity, the value of 6 in (11) 
can be made arbitrarily small based on the choice of e. This in turn imphcs that 
the two eigenvalues can be made arbitrarily close to each other with e. On the 
other hand, for a fixed M, as e tends to zero, (11) becomes 



4.2 Frames in C^; d > 2 



For general d and Af , in order to use existing results on the conccintration of 
eigenvalues of random matrices [23, 24], a slightly different construction of the 
frame needs to be considered. Let {Ymn]m,n& be i.i.d. random variables 
following a Gaussian distribution with mean zero and variance . It can be 
shown that 

E{e « ' "") = e 2 U ; 

and the variance 

One can define the following two dimensional sequence. For m,n gZ, 

Xmn = e « - e 2 U ; . 
Consider the mapping u : Z — )• given by 

/ Xu \ 



X2l 



V Xm j 



(13) 



As before, let M > cZ and consider the set of M unit vectors V = {v{l), v{2), . . . , v{M)} 
in C*. The frame operator of this set is 



Xu Xi2 

X21 X22 



Xdl Xd2 



XiM 

X2M 



Xd_ 



M 



Xu X21 

X\2 X22 



Xdl 
Xdl 



X 



IM X2M 



X. 



dM 
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Let 



Vd 



Xll Xi2 

X21 X22 



Xdl Xd2 



XiM 

X2M 



XdM 



(14) 



so that = AA* .^The matrix A has entries with mean zero and variance 
(7^ = A(l-e~'^^(^) ). According to results in [23], if ^g- -> c as d, M 00, then 



the smallest and largest eigenvalues of converge almost surely to a-^(l • 
and (7^(1 + y/c)"^ respectively. 



Theorem 4.5. Let si{A) < S2{A) < . . . < sd{A) be the singular values of the 
matrix A given by (14)- Then the following hold, 
(a) Given eo, there is a large enough d such that 



P Sd{A) > a 1 




P{si{A)<c^)<e-^'^ 
where c\ and C2 are universal positive constants. 



(15) 



(16) 



Proof. Let Sd be the mapping that associates to a matrix A it largest singular 
value. Equip C^*^ with the Frobenius norm 

\\Af := Tr{AA*) = J2\^mnf- 



Then the mapping Sd is convex and 1-Lipschitz in the sense that 

\sd{A) - SdiA')\ < \\A-A'\\ 

for all pairs {A, A') of d by M matrices [24]. 

We think of A as a random vector in R^"^*^ . The real and imaginary parts of 
the entries of are supported in [— -^j- Let P be a product measure on 
[— -i=, -i=]2dM Tjjgjj g^g consequence of concentration inequality (Corollary 4.10, 
[24]) we have 

P{\sd{A)-m{sd)\>r)<Ae-'-"''^'^ 

where m{sd) is the median of Sd{A). It is known that the minimum and maxi- 
mum singular values of A converge almost surely to — \/c) and (t(1 + ^/c) 
respectively as rf, M tend to infinity and -p^ — )■ c. As a consequence, for each 
Co and M sufficiently large, one can show that the medians belong to the fixed 

interval [a{l - J jj) - eo, <7(1 + V ig) + eo] which gives 




< 2e-'-'''/i^ 
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For the smallest singular value we cannot use the coneentration inequality as 
used for Sd since the smallest singular value is not convex. However, following 
results in [25] (Theorem 3.1) that have been used in [26] in a similar situation 
as here, one can say that whenever M > {l + 5)d where 6 is greater than a small 
constant, 

P{si{A) < ci) < e-^'^ 
where ci and C2 are positive universal constants. □ 

Remark 4.6. Note that the square of the singular values of A are the eigen- 
values of T and so the estimates given in (15)-(16) give insight into the corre- 
sponding deviation of the eigenvalues of the frame operator T. 

Remark 4.7 {Connection to compressed sensing). The theory of compressed 
sensing [27, 28, 29] states that it is possible to recover a sparse signal from a small 
number of measurements. A signal x G C*'^ is fc-sparse in a basis ^ = {tpj}j^i 
if X is a weighted superposition of at most k elements of Compressed sensing 
broadly refers to the inverse problem of reconstructing such a signal x from 
linear measurements {yi — {x^(j)i) \ (. = !,...,£?} with d < M, ideally with 
d M. In the general setting, one has ^*x = y, where $ is a d x M sensing 
matrix having the measurement vectors as its columns, a; is a length-M signal 
and y is a length- d measurement. 

The standard compressed sensing technique guarantees exact recovery of 
the original signal with very high probability if the sensing matrix satisfies the 
Restricted Isometry Property (RIP). This means that for a fixed k, there exists 
a small number 6k, such that 

{l-Sk)Ml<\\<Px\\l<{l + 5k)\\x\\l, 

for any fc-sparse signal x. By imitating the work done in [26] (Lemmas 4.1 and 
4.2), it can be shown, due to Theorem 4.5, that matrices A of the type given 
in (14) satisfy the RIP condition and can therefore be used as measurement 
matrices in compressed sensing. These matrices are different from the traditional 
random matrices used in compressed sensing in that their entries are complex- 
valued unimodular instead of being real-valued and not unimodular. 

Example 4.8. This example illustrates the ideas in this subsection. First 
consider M = 5 and d = 3 so that there are 5 vectors in C'^. Taking from a 
normal distribution with mean and variance cr = 1, a realization of the matrix 

[^Ti]l<m<3,l<n<5 is 

" -0.0353 0.5004 -0.6299 -0.1472 0.4003 " 
-0.4804 -0.9344 0.4220 -0.9509 0.2783 . 
-0.8609 -0.4822 -0.4680 -0.0913 1.2284 

Then taking e = 0.001, A = -Ue^'^-^^] is 



V3 



-0.27- 0.96i 
-0.92 - 0.39i 
0.74 -I- 0.68i 



-0.89 + 0.44i 
-0.89 - 0.46i 
0.16-0.99i 



0.85 + 0.52i 
0.99 -I- 0.05i 
0.99 -I- 0.09i 



0.33 - 0.94i 
0.93 -I- 0.37i 
-0.30 - 0.95i 



-0.24+0.97i 
-0.47-F0.88i 
-0.74-F0.67i 
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Figure 1: Behavior of the condition number of the frame operator with increas- 
ing size of the frame; e — 0.0001, d=2>,a — \ 



The condition number, ratio of the maximum and minimum eigenvalues, of 
T — 4.8667. As the number of vectors M is increased, the condition number 
gets closer to 1. Figure 1 shows the behavior of the condition number with the 
increase in the number of vectors. 



5 Conclusions 

The construction of discrete unimodular stochastic waveforms with arbitrarily 
small expected autocorrelation has been proposed. This is motivated by the 
usefulness of such waveforms in the areas of radar and communications. The 
family of random variables that can be used for this purpose has been character- 
ized. Such construction been done in one dimension and generalized to higher 
dimensions. Further, such waveforms have been used to construct frames in 
and the frame properties of such frames have been studied. Using Brownian 
motion, this idea is also extended to the construction of continuous unimodular 
stochastic waveforms whose autocorrelation can be made arbitrarily small in 
expectation. 
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